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[57] ABSTRACT 

A set of parameters (W) is calculated for an equalizer which 
equalizes a transmission line impulse response (UR) in such 
a way that the equalized impulse response (EIR) approxi- 
mates a target impulse response (TIR) with predetermined 
length (L). An error function (E) is minimized via 
eigenvalue -eigenvector calculation of a channel dependent 
matrix (O). This channel dependent matrix (O) is built up 
from a signal (X), a distorted signal (Y), and length (L) and 
delay (D) of the target impulse response (TIR). The error 
function (E) comprises a first contribution which represents 
the difference between the equalized impulse response (EIR) 
and the target impulse response (TIR) and a second contri- 
bution which represents the energy (E^) that is transmitted 
in unused frequency bands. The eigenvector corresponding 
to the minimum eigenvalue of the channel dependent matrix 
(O) represents the set of equalizer parameters (W). 

9 Claims, 5 Drawing Sheets 
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METHOD FOR TRANSMISSION LINE 
IMPULSE RESPONSE EQUALIZATION AND 
A DEVICE TO PERFORM THIS METHOD 

TECHNICAL FIELD 

The present invention relates to a method for calculating 
a set of equalizer parameters, a device for calculating such 
a set of equalizer parameters, and a demodulator including 
such a device. 

BACKGROUND OF THE INVENTION 

Such a method and equipment to perform such a method 
are already known in the art, e.g. from the article ' Adaptive 
Channel Memory Truncation for Maximum Likelihood 
Sequence Estimation', written by D. D. Falconer and R R. 
Magee, Jr., and published in The Bell System technical 
Journal, Vol 52, No. 9, November 1973, from page 1541 to 
page 1562. As described in the second paragraph of the 
introduction on page 1541, an equalizer or linear prefilter 
which forces the overall impulse response of the channel 
prefilter combination to approximate a desired truncated 
impulse response is presented therein. The overall impulse 
response of the channel prefilter combination is equivalent 
to the equalized impulse response mentioned in the pre- 
ambles of the independent claims hereof, while therefor the 
desired truncated impulse response in the above mentioned 
article is called a target impulse response with predeter- 
mined length in the present application. In the known 
method, the desired impulse response is obtained by calcu- 
lating an eigenvector corresponding to the minimum eigen- 
value of a certain channel dependent matrix. This channel 
dependent matrix and more specifically the minimum eigen- 
value thereof is representative for a mean square error 
between the output of the prefilter and the output of the 
desired prefilter, the desired prefilter having an impulse 
response equal to the just mentioned desired truncated 
impulse response. When minimizing the mean square error, 
the desired impulse response is optimized. As is proven in 
section III of the above cited article, the known method then 
adaptively or iteratively calculates the linear prefilter param- 
eters from the desired truncated impulse response. Thereto, 
in successive iterations, the Viterbi algorithm is applied to 
the prefilter output, the result thereof is applied to the desired 
prefilter and the output of this desired prefilter is subtracted 
from the prefilter output in such a way that an error sequence 
for adjustment of the desired impulse response and prefilter 
parameter set is obtained. The receiver structure allowing to 
perform the known method is drawn in FIG. 3 on page 1545 
of the above mentioned article. Although the known tech- 
nique or derivatives thereof wherein e.g. the Viterbi algo- 
rithm is not used exhibit very low mean square error values 
and as a consequence a theoretically negligible intersymbol 
interference, its applicability in ADSL (Asymmetric Digital 
Subscriber Line) is doubtful. The main reason is linked to 
the finite precision limitations of the hardware. Indeed, the 
equalizer transfer function obtained when applying the 
known technique in ADSL; will boost the unused frequency 
band, i.e. the frequency band below 25 kHz, and will 
attenuate the tones in the ADSL passband, which approxi- 
mately cover the frequency band from 25 kHz to 138 kHz. 
A first consequence is the decrease of numerical stability. It 
is even possible that the attenuation is so severe that the first 
stages of an FFT demodulator as is used in ADSL receivers 
rounds off the frequency components to zero such that the 
frequency domain equalizer of the receiver is unable to 
boost these tones again. An additional consequence is that, 
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by implementing the known method in ADSL, the attenua- 
tion is very severe at the filter boundaries between the 
unused and used frequency bands. This is a result of the very 
sharp filter edges there. The pilot tone will very likely be 
located close to the region where the transmission channel 
exhibits a relatively low attenuation. This is mostly located 
near the lower downstream filter edge. Synchronization may 
therefore be subjected to an unexpected decrease of the 
signal to noise ratio of the pilot tone. 

DISCLOSURE OF INVENTION 

An object of the present invention is to provide a method 
and equipment for calculating a set of equalizer parameters 
of the above known type but wherein numerical instability 
due to amplification of unused frequency bands and attenu- 
ation of used frequency bands is optimized, and wherein the 
above mentioned synchronization problem is avoided. 

According to the invention, this object is achieved in the 
method for calculating a set of T equalizer parameters W, T 
being a positive integer number, for an equalizer included in 
a receiver to which a transmitter is coupled via a transmis- 
sion line and which is able to receive a distorted signal Y 
distorted due to transmission over said transmission line of 
a signal X, and which is adapted to equalize a transmission 
line impulse response LIR of said transmission line and to 
thereby generate an equalized impulse response EIR which 
approximates a target impulse response TIR with a prede- 
termined target impulse response length L delayed over a 
predetermined delay D, said method comprising the steps of: 

a. transmitting said signal X and determining said dis- 
torted signal Y in said receiver; 

b. defining a first integer number of sample periods which 
represents said target impulse response length L; 

c. defining a second integer number of sample periods 
which represents said predetermined delay D; 

d. generating a channel dependent matrix O from said 
signal X, said distorted signal Y, said target impulse 
response length L and said predetermined delay D, said 
channel dependent matrix O being related to an error 
function that has to be minimized; 

e. calculating a minimum eigenvalue K MJN o^ said channel 
dependent matrix O; 

f. calculating an eigenvector of said channel dependent 
matrix O, said eigenvector corresponding to said mini- 
mal eigenvalue Xj^; 

g. generating said set of T equalizer parameters W from 
said eigenvector, characterized in that said error func- 
tion comprises a first contribution, representing a dif- 
ference between said equalized impulse response EIR 
and said target impulse response TIR, and a second 
contribution representing an amount of energy E DC 
transmitted in unused frequency bands. 

According to the invention, the above object is also 
achieved in a device for calculating a set of T equalizer 
parameters W, T being a positive integer number, for an 
equalizer included in a receiver to which a transmitter is 
coupled via a transmission line and which is able to receive 
a distorted signal Y distorted due to transmission over said 
transmission line of a signal X, and which is adapted to 
equalize a transmission line impulse response LIR of said 
transmission line and to thereby generate an equalized 
impulse response EIR which approximates a target impulse 
response TIR with a predetermined target impulse response 
length, when delayed over a predetermined delay D, said 
device including: 
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a. a determining means to determine said distorted signal c. a second definition means to define a second integer 

Y received upon transmission of said signal X; number of sample periods which represents said pre- 

b. a first definition means to define a first integer number determined delay D; 

of sample periods which represents said target impulse d. a first generation means to generate a channel depen- 

response length L; 5 dent matrix O from said signal X, said distorted signal 

c. a second definition means to define a second integer X said target impulse response length L and said 
number of sample periods which represents said pre- predetermined delay D, said channel dependent matrix 
determined delay D; 0 being related to an error function which has to be 

d. a first generation means to generate a channel depen- minimized; 

dent matrix O from signal X, said distorted signal Y, 3 o e - a firet calculation means to calculate a minimal eigen- 

said target impulse response length L and said prede- value of said channel dependent matrix 0; 

termined delay D, said channel dependent matrix O f , a second calculation means to calculate an eigenvector 

being related to an error function which has to be 0 f sa id channel dependent matrix O, said eigenvector 

minimized; corresponding to said minimal eigenvalue ~^ M in> arj d 

e. a first calculation means to calculate a minimal eigen- 15 g a second generation means to generate said set of T 
value \ MIN of said channel dependent matrix 0; equalizer parameters W from said eigenvector, charac- 

f. a second calculation means to calculate an eigenvector terized in that said error function comprises a first 
of said channel dependent matrix O, said eigenvector contribution, representing a difference between said 
corresponding to said minimal eigenvalue X MIN ; and equalized impulse response EIR and said target impulse 

g. a second generation means to generate said set of T 2 o response TIR, and a second contribution representing 
equalizer parameters W from said eigenvector, charac- an amount of energy £ DC transmitted in unused fre- 
terized in that said error function comprises a first quency bands. 

contribution, representing a difference between said Indeed, by adding to the difference between the equalized 

equalized impulse response EIR and said target impulse and the target impulse response, a second contribution, an 

response TIR, and a second contribution representing 2 s additional degree of freedom is provided which allows 

an amount of energy E DC transmitted in unused fre- shaping of the equalizer impulse response. By adding a new 

quency bands. error source to the mean square error criterion, the present 

According to the invention, the above object is also invention prevents boosting unused frequency bands instead 

achieved in a demodulator DEMOD including between an 0 f used frequency bands. The higher the energy located in 

input DI and an output DO thereof, the cascade connection 30 the unused frequency band, the higher the additional error, 

of an analog to digital converter ADC, a time domain As a result, the present method searches for equalizer 

equalizer TEQ, a serial to parallel converter SPC a fast parameters which define an equalizer impulse response that 

fourier transform processing unit FFT, a frequency domain does not boost unused frequency bands. From the above 

equalizer FEQ and a demapper DM, said analog to digital ' explanation, it is clear that in the present invention, the 

converter ADC being provided to convert a disturbed input 35 shortcomings with respect to numerical stability and syn- 

signal applied to said demodulator input DI into a time chronization are overcome. 

domain digital signal, said time domain equalizer TEQ being a further feature of the present invention is that the first 

provided to digitally filter said time domain digital signal contribution to the channel dependent matrix, representing 

and to thereby generate a time domain equalized signal, said the difference between equalized and target impulse 

serial to parallel converter SPC being provided to transform 40 response, has the form R M ^-R Kr RjKy" 1 .R K /- 

said time domain equalized signal into a time domain Indeed, as will be proven later on in the description, 

parallel sequence, said fast fourier transform processing unit modelling the transmission line in a certain predetermined 

FFT being provided to process said time domain parallel wav , allows the first part of the channel dependent matrix 0 

sequence and to thereby generate a frequency domain par- as given by the to build up expression to be denied, 

allel sequence, said frequency domain equalizer being pro- 45 [ n different implementations of the present method, step 

vided to digitally filter said frequency domain parallel d of generating a channel dependent matrix O may be 

sequence, and, said demapper DM being provided to simplified. Therefore, assumptions can be made. Assuming 

demodulate said frequency domain parallel sequence and to that the transmission line is not plagued by noise, or 

generate therefrom an output signal applied to said demodu- assuming, that the spectral density of the input signal applied 

lator output DO, said time domain equalizer TEQ including 50 to the transmission line by the transmitter has a flat shape, 

a device for calculating a set of T equalizer parameters W, permits amplification of the mathematical expression of the 

T being a positive integer number, for said time domain channel dependent matrix O, and reduces the number of 

equalizer TEQ, said demodulator DEMOD being coupled to calculations necessary in step d to build up this matrix. For 

a transmitter via a transmission line and being able to receive further details, reference is made to the description later on 

a distorted signal Y distorted due to transmission over said 55 in this specification. 

transmission line of a signal X, and which is adapted to An additional characteristic feature of the present inven- 
equalize a transmission fine impulse response LIR of said uon i s that, in implementations thereof wherein the trans- 
transmission line and to thereby generate an equalized mission line is not supposed to be noiseless, the noise is 
impulse response EIR which approximates a target impulse measured in a step included in the present method before the 
response TIR with a predetermined target impulse response 60 channel dependent matrix O is built up. 
length L> when delayed over a predetermined delay D, said a further characteristic feature of the present method is 
device including: that, in existing implementations thereof wherein the trans- 

a. A determining means to determine said distorted signal mission fine is not supposed to be noiseless, the noise is 

Y received upon transmission of said signal X; taken into account. Indeed, by modelling noise measured on 

b. a first definition means to define a first integer number 65 the transmission fine via an additive noise source, connected 
of sample periods which represents said target impulse in cascade to a noiseless transmission line, elements of the 
response length; autocorrelation matrix RYY are modified as described 
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above. Later on, it will be seen how the expression of RYY pared to the length of data symbols transmitted over this 

is mathematically derived. transmission line. In known multicarrier systems, however, 

Yet a further characteristic feature of the present method acceptable latency or maximum acceptable delay between 
is that, the predetermined delay (D) depends upon the input and output of data symbols transmitted over a trans- 
amplitude and the index number of a sample of said trans- 5 mission line limits the length of these symbols. As is 
mission line impulse response (LIR) with maximum ampli- specified in section 6.2.11 on page 22 of the draft American 
tude. The second number of sample periods can be National Standard for Telecommunications on ADSL 
expressed as a function of the amplitude and the sample (Asymmetric Digital Subscriber line) published b> ' ANSI 
index number of the transmission Hue impulse response ^rnencan National Starulards tetitute) m Ajpnl 1994 the 
sample with the highest amplitude. This feature corresponds 10 ^ of T tone ^ M T> s 5™bol * bmited. In 
to the intuitive feeling that the shorter target impulse ^ applications, the transmission hue impulse response 
response will approximate the transmission line impulse len S th * not negligible compared to the symbol length and 
response in an optimized way only if this shorter target *tersymbol interference wnl occur. As a consequence suc- 
impulse response is delayed in such a way that the location cessiv f ^ f nsmitted symbols are spread out and interfere as 
in time of its peak values approximate the peak values of the 15 a result of whlch the transmission is distorted, 
transmission line impulse response. ^ general, two different techniques are known to com- 

In a further implementation of the method according to pensate for such intersymbol interference. In a first 
the present invention, before executing steps f and g, steps technique, called cyclic prefix extension, data symbols are 
c, d and e are repeatedly executed for a sequence of extended at the transmitters. Gaps are created at the edges of 
acceptable delay values. The second integer number of 20 the data symbols e.g. by inserting spare time between 
sample periods is different in successive executions of step successive symbols. In these gaps intersymbol interference 
c but laying within an interval of acceptable delay values, can occlir without distorting the transmission. The transmis- 
said channel dependent matrix O generated in successive sion efficiency or throughput efficiency, however, decreases 
executions of step d as a result thereof being different in such increasing lengths of the gaps or cyclic extensions. In 
a way that a plurality of minimum eigenvalues similar to 25 lhe sec °nd intersymbol interference compensating 
said minimum eigenvalue \ MIN is calculated in successive technique, the transmission line impulse response is equal- 
executions of said step e, and that to perform steps f and g, b ? m equalizer, which is a digital filter at the receiver's 
only a said channel dependent matrix O with smallest said Sldc - ^ objective of the equalizer is to compress the 
minimum eigenvalue X uar is selected. In this way, for each transmission line impulse response to a shorter overall 
delay value within an interval of acceptable delay values, a 30 impulse response. Full equalization of the transmission line 
corresponding channel dependent matrix having its specific ^pulse response however requires a large amount of corn- 
minimum eigenvalue is found. It will again be proven later pntation. The draft ADSL Standard already mentioned above 
on that this minimum eigenvalue is a measure for the error specifies the combination of both techniques. Equalizing the 
function that must be minimized, and therefore, in this transmission line impulse response renders the transmission 
further embodiment, only the channel dependent matrix with 35 llDe impulse response length short enough to enable 
the smallest minimum eigenvalue is maintained to calculate intersymbol interference compensation via cyclic prefix 
therefrom in steps f and g the set of equalizer parameters. extension with tolerable prefix length, and vice versa: com- 

These and other objects, features and advantages of the pensation via cyclic prefix extension allows the equalizer to 

present invention will become more apparent in light of the be of acceptable computation complexity, 

following detailed description of a best mode embodiment 40 In the following paragraphs, an algorithm will be derived 

thereof, as illustrated in the accompanying drawings. wm *ch enables calculation of equalizer parameters in a fast 

and stable way. Before entering the mathematical approach 

BRIEF DESCRIPTION OF THE DRAWINGS of the problem, the ADSL (Asymmetric Digital Subscriber 

Line) environment, i.e. the demodulating part, wherein such 

FIG. 1 represents a block scheme of an embodiment of a 45 m equalizer is used, will be described, 

demodulator according to the present invention; and The demodulator DEMOD 10 of FIG. 1 includes an 

FIG. 2 represents a digital model for the cascade connec- analog to digital converter ADC 14, a time domain equalizer 

tion of a transmission line and an embodiment of an equal- TEQ 16, a serial to parallel converter SPC 18, a fast fourier 

izer whose time domain equalizer parameters are calculated transform processing unit FFT 20, a frequency domain 

according to the present method. 50 equalizer FEQ 22, and a demapper DM 24 which are 

FIGS. 3 A and 3B represent a flowchart for the method of connected in cascade between an input DI 26 and an output 

calculating an set of equalizer parameters (W). DO 28 of the demodulator DEMOD 10. 

FIG. 4 represents a block diagram of a device for calcu- Distorted signals at the output of the transmission line, 

lating a set of equalizer parameters (W). connected to the demodulator input DI 26 but not shown in 

FIG. 5 represents a block diagram of the time domain 55 eater the demodulator DEMOD 10 via the input DI 

equalizer TEQ ^ ucn a distorted signal is sampled by the analog to digital 

„ r ' T , , , converter ADC 14 to be converted into a time domain digital 

The reference numeral for a particular step of the method , . ,. r ? 

i_ • ™^ ia » prrn ,r.u u ~i • i sequence which is applied to the equalizer TEQ 16. The time 

shown in FIG. 3A and FIG. 3B has been made equivalent to , _ . «. „ ™~ 1<: . n ?• t ni . „ r , _ . nf 

e . - „ , ? i . • . domam equalizer TEQ 16 is a digital filter whose set or 

the reference number of the corresponding block which 4 . i i ♦ j u a t. ij j * ■ f 

. . , . , r . A 1 60 parameters is calculated based on channel dependent lnior- 

performs that step m the device shown in FIG. 4. ^ accordance with ^ method of the present mven . 

BEST MODE FOR CARRYING OUT THE tion * In by ^ally mterm S the time domain digital 

INVENTION sequence, the overall impulse response length is shortened. 

The resulting equalized time domain digital sequence is 
In a theoretically ideal muliticarrier system, i.e. in a 65 transformed into an equalized parallel time domain-digital 
system with infinitely large data symbols, the impulse _sequence.by_the^serial to-pa^aUel-converte^SPG-18.b^ last " i 
response length of the transmission line is negligible com- \ fourier transform processing unit FFT 20jx>nverts the^uaF ' / 
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^Ized parallel time domain digital sequence into an equalized 
frequency domain parallel digital sequence and applies this 
sequence to the frequency domain equalizer FEQ 22 which 
restores amplitudeiM phisfc of the transmitted vectors, and 
tojhe demappjer DM 24.-This equalized^frequency domain 
parallel digital sequence is a sequence of real data which 
represents a set of modulated carriers. These carriers have") 
equidistant frequencies, carry different amounts of data bits 
which are modulated thereon via different modulation tech- 
niques. To determine the number of data bits that are . 
allocated to different carriers and to determine the modula- 
tion techniques used to modulate different carriers, a mapper 
included in the transmitter must communicate with the 
demapper DM 24 in the receiver. This mapper e.g. assigns 
2 data bits to the first carrier and decides to modulate these 
2 bits thereon via 4 QAM modulation, assigns 4 data bits to : 
the second carrier and decides that these 4 data bits are 
modulated thereon via 16 QAM modulation, and so on. The 
demapper DM 24 included in the demodulator DEMOD 10 
at the receiver's side is made aware of the modulation 
techniques used for the different carriers and the number of 
bits allocated to each of the carriers. The task of the 
demapper DM 24 is thus to demodulate each carrier and to 
obtain therefrom the exact number of data bits modulated 
thereon. The demapper DM 24 then constitutes an output 
sequence and supplies this output sequence to the demodu- 
lator output DO 28. 

The block scheme of the Discrete Multi Tone demodulator 
DEMOD 10 in FIG. 1 will not be described in further detail 
since such a detailed description is of no importance for the 
present invention. Further details with respect to ADSL 
requirements are described in the already mentioned draft 
ANSI Standard on ADSL while specific implementations of 
Discrete Multi Tone modulators and demodulators are found 
in the articles ' Amulticarrier El -HDSL Transceiver System 
with Coded Modulation' , written by Peter S. Chow, Naofal 
Al-Dhahir, John M. CiofB and John A. C. Bingham and 
published in the issue Nr. 3 May/June 1993 of the Journal 
of European Transactions on Telecommunications and 
Related Technologies (ETC), pages 257-266, and Perfor- 
mance Evaluation of a Multichannel Transceiver System for 
ADSL and VHDSL Services' from Peter S. Chow et al., 
published in the issue Nr. 6 August 1991 of the Journal of 
European Transactions on Telecommunications and Related 
Technologies (ETT), pages 909-919. 

In the following paragraphs, the method according to the 
present invention will be mathematically derived. Starting 
from a time domain model, a mathematical approach of the 
problem will lead to an expression of an error function and 
a channel dependent matrix which form the basis of the 
present method. 

The time domain model is drawn in FIG. 2. Therein, the 
transmission line is represented by the cascade connection of 
a noiseless transmission line with impulse response LIR 60 
and an additive noise source wherein noise N 76 is added via 
a summer SM 62 to the output signal of the noiseless 
transmission line II 92. This output signal enters the summer 
SM 62 via a first input II 92, while the noise N 76 enters the 
summer SM 62 via a second input 1290 thereof. An undis- 
torted signal X 77 applied to the noise N 76 transmission line 
appears at the output thereof as a distorted signal Y 78, 
distorted due to the transmission line impulse response LIR 
60 and due to the presence of the noise N 76. To compensate 
for these disturbances, a digital filter 64 filters the distorted 
output signal Y 78. This digital filter 64 is included in an 
equalizer 90 characterized by a set of parameters W. The 
transmission line impulse response LIR 60 and the equalizer 
parameters W constitute an equalized impulse response EIR 
72. 



20 



As already mentioned when describing the problem of 
intersymbol interference above the equalized impulse 
response EIR 72, should be short in length to enable 
avoiding intersymbol interference. Therefore, the equalized 
impulse response EIR 72 in the digital model of FIG. 2 is 
compared to a target impulse response TIR 68, not further 
specified, but with a predetermined length L representing the 
guardband length or cyclic prefix length. This target impulse 
response TIR 68, delayed over a delay D 104 must approxi- 
mate the equalized impulse response EIR 72 to guarantee 
negligible intersymbol interference. Consequently, the mis- 
match between the equalized impulse response EIR 72 and 
the target impulse response TIR 68 must be minimized. As 
will be seen in the calculations later on, an error E 88 
between the output 71 80 of the equalizer 90 and a fictional 
output 72 84 of the target impulse response TTR 68 consti- 
tutes the first contribution to an error E 88 that must be 
minimized. These outputs, Zl 80 and 72 84 in FIG. 2 are 
applied to a non-inverting input NI 82 and an inverting input 
II 86 of a subtracter SB 70 respectively. 

The following definitions are used in the mathematical 
approach in the following description: 



25 



30 



35 



2 



45 



50 



55 



60 



65 



Xh,k" 



Xb,k m 



Y„,k- 



hi 



b 0 
bi 

bL-i 



X k . t 

Xk-D 
Xjt-D-l 

Xk-D-L+V 

n 
y*-i 



the Z- trans form A (2) of a digital sequence a,- 
whose index » covers a range from -oo to +00; 

the length of the transmission line 

impulse response LIR; 
the number of taps in the digital filter of the 
equalizer or, in other words, the number of 
time domain equalizer parameters W\ 

the length of the target impulse response 
TIR; 

the delay applied to the target impulse 
response TIR; 



the line impulse response vector; 



the equalizer time domain vector; 



the target impulse response vector 



a signal vector related to the line impulse 
response; 



a signal vector related to the target impulse 
response; 



a distorted signal vector related to the 
equalizer time domain; 
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-continued 

Wopt : a a optimal equalizer time domain vector; 
Bopt : an optimal target impulse response vector; 
and 

O : a channel dependent matrix. 

As shown in FIG. 2 when an undistorted signal X 77 is 
applied to the transmission line, a distorted signal Y 78 
appears at the output thereof. The k'th sample, k being an 
integer, of the distorted signal Y 78 does not only depend 
upon the k'th sample of the undistorted signal X 77 but also 
upon M-l previous samples of this undistorted signal X 77. 
This is due to the transmission line impulse response LIR 60 
which has length M. The transmission line is not perfect and 
as a result input samples transmitted in the past still have ao 
influence on future output samples. The k'th output sample 
Y* of the transmission line can be expressed as follows: 



- E h r x k -i+n k 
1=0 



(1) 



Herein, i represents an index and the influence of noise od 
the transmission line is concentrated in the term n^ this 
being a result of the digital model of FIG. 2. Using the above 
denned vectors, the expression in (1) can be simplified as 
follows: 



y k -H r X Kk +n k 



(2) 



Herein, H r represents the transposed vector of H. 

The distorted output signal Y 78 is applied to the equalizer 
90 to be digitally filtered thereby. Since the digital filter 64 
has T taps, the k'th sample zl* of the equalized output signal 
Zl 80 is given by the convolution 



2\k> 



r-i 

> 2 wryk-i 
1=0 



or, when using vector expression: 



(3) 



(4) 



35 



L-l 
1=0 



(5) 



10 



Therein, B r denotes the transposed vector of B. 

As shown in FIG. 2 the target impulse response TIR 68, 
the time domain equalizer parameters W and the delay D 
104 now have to be determined in such a way that the 

5 equalized impulse response EIR 72 approaches the target 
impulse response TIR 68 as closely as possible. Therefore, 
output signals Zl 80 and Z2 84 are compared and the 
mismatch between both is minimized. This error E 88 is 
mathematically expressed by a sequence E whose k'th 

10 component c k is given by: 



/=0 



L-l 

2 brx k -D-i 
1=0 



CO 



15 or, when using the vector notation: 

e k -w T Y Wik -B T x bJ , (8) 
Thus, the error function that has to be minimized is equal to: 



20 



E{e?}-E{^Y wjr B T X bJ f} 



(9) 



Therein, E 88 denotes the mean square error between Zl 80 
and Z2 84. This mean square error is not dependent upon k 
and can thus be rewritten as follows: 



25 



E{e 2 }=E{(W r Y tt/ -B r X b ) 2 } 



(10) 



Since 



W r represents the transposed vector of W herein. 

As shown in FIG. 2 when the undistorted signal X 77 is 
applied to a fictional filter or fictional transmission line with 
a target impulse response TIR 68 whose length is restricted 45 
to L taps, a second, fictional, output signal Z2 84 is obtained. 
The target impulse response TIR 68 is defined by the target 
impulse response vector B. The shape of the equalized 
impulse response EIR 72, which comprises the line impulse 
response LIR 60 and the digital filter in the time domain 50 
equalizer TEQ 16 of FIG. 1, response will be determined in 
such a way that it approaches the target impulse response 
TIR 68. To allow the equalized impulse response EIR 72 to 
approach the target impulse response TIR 68, a delay D 104 
is applied to the target impulse response TIR 68. This delay ss 
D 104 should be considered as an additional parameter 
which must be determined somewhere in the method as 
shown in step 104 of FIG. 3A and block 104 of FIG. 4. The 
delayed target impulse response transforms the undistorted 
signal X 77 into the second output signal Z2 84, whose 60 
sample with index k is given by the convolution: 



(W r .Y„-B T X b ) 2 =W r .Y tv .Yj.W-ZW T . Y„X b T £+B T JC b X b T £ (11) 

and by introducing the following correlation matrices: 

Ryy-E{Y„Yj) 

Rrx-EiY^r} (12) 

(10) can be rewritten in the form: 

E{e 2 }-W r Jiyy.W+B T Jl Jtx £-2B T R YX r .W (13) 

With the following substitutions 



Rxx= 



Ryy- 



Ryx- 



RrM /Ui) ... lUL-i) 

*~(D 

R„(L - 1) ... Rjl) *„(0) 

J?„(D) R„{1) ... Ryy(T~l) 

*>y(°) 



*»<r-i) ... 

Ry:(D-T+l) 



Ryx(P + 1) 

M0) 



(14) 



*>y(U) 

.. RyxiD+L-l) 

.. Ryi(p + L-T) 



If it is assumed that the noise N 76 is uncorrelated with the 
undistorted signal X 77 then the following results from (14). 



or in matrix notation: 



65 



7H k "B T Ji bwk 



(<0 



(15) 
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12 



It is remarked here that: 



(m-i \ 

= E< X hi Xi+t.i xi V 

I " J 

2 hr E{xuk-i ■ 
M) 

M-l 

= R„(k)*h k 
and similar calculations prove that: 

Ry^JkYihfhd (17) 

To minimize the mean square error of (13) zero crossings 
of a partial differential should be found. The target impulse 
response vector B which permits the minimization of the 
mean square error, i.e. the optimal target impulse response 
vector Bopt, in (13) is obtained by solving the following 
equation: 



O^Ryy-R yx^xx' 1 -Ryx T (25) 

^ When minimizing E{e 2 } by searching the equalizer time 
domain vector W which minimizes W r .O.W it is obvious 
_ that the solution: 



10 



15 





0 




0 




0 


and thus: Bopt - 


0 




0 




0 



(26) 



will be found. However, this solution would imply that no 
signal could be transmitted over the transmission line. To 
obtain another solution, different from the zero -filter, an 
energy constraint: 



(27) 



has to be fulfilled. 
Thus, the mean square error (24) must be minimized by 
taking into account the constraint (27). This is done by 
20 solving the equations: 



6(B) 



(18) 



25 



fi(E{«2}-A.-(WT. w- 1)) 

bw 

fi(g{c a }-A,-(W r -W-l)) 



-0 



= 0 



(28) 



Substitution of (13) in (18) results in the following: 

2R XX £-2R VX .W-Q (19) 



Substituting (24) in (28) and calculating the partial deriva- 
tive results in: 



2- o- w-2-\- w=o 
m- w-i oo 



(29) 



30 

which implies that the optimal target impulse response 
vector Bopt can be calculated once the equalizer time 
domain vector W is known via the formula: 

Hie second equation in (29) again represents constraint (27) 
Boptwtjfl^-Anr ( 20 ) 35 which must be fulfilled. The first equation is equivalent to: 

Herein, Rxy -1 represents the inverse matrix of the above 

defined autocorrelation matrix R^ and Ryx T represents the o.w*\.w (30) 

transposed matrix of the crosscorrelation matrix R^. 

To obtain an expression of the mean square error E which The equation (30) is an eigenvalue-eigenvector problem, 
depends only upon the equalizer time domain vector W, (20) 40 Since 
is substituted in (13). Thus: 

min(£{e 2 }) - minfW 7 ' ■ O ■ W) - min(^ r • k - W) = (31) 

E{<?)=W J JlyrW+W 7 . R^ {K^-yjtjarR^Jt^.W- 

2W r Jt YX .Qixx-YAYx r -W (21) ' W ^ * mm W> 



which is further simplified as follows: 

E{<?} « WiRYY + RYxiRixY'RYx- (22) 

2Ryk(R x S ( ) t 'RI x )W 

- Wprr-RniR-AVRidW 

From Rxx-Rxx~ 1= ^> wherein 1 represents the unity matrix, it 
follows that (Kxx'Y^x/^h and thus that (R JHr - 1 ) r - 

As a result, (22) can be transformed into: 



45 



50 



55 



E{e?}=W r .(R Y y-R YX .<R xx 7 )~RYX T )W 



(23) 



From the definition of the autocorrelation matrix in 6Q 
(14), it further follows that R^^Rjer ^ a result of which: 



E{c 2 } = ^(RrY-RYxR-A-R^W 

- w? o w 



(24) 



Herein, O is a channel dependent matrix given by the 
expression: 



the optimal equalizer time domain vector Wopt equals the 
eigenvector of O which corresponds to the minimum eigen- 
value Xmin thereof. 

Thus, summarizing the above calculations as shown in 
FIG. 4, the channel dependent matrix O must be built up as 
shown in block 130, the minimum eigenvalue Xmin as 
shown in block 132 and corresponding eigenvector of 0 134 
must be calculated. This eigenvector equals the optimized 
equalizer time domain vector Wopt and finally, the optimal 
target impulse response vector Bopt is calculated from Wopt 
via equation (20) as shown in block 138. Once the matrix O 
is known, all steps to derive therefrom the equalizer time 
domain parameters W are thus known. 

Formula (25) shows that the channel dependent matrix O 
only depends on the undistorted signal X 77 and the received 
distorted signal Y 78. The undistorted signal X 77 is trans- 
mitted upon initialization of the system and is well known by 
the receiver. At reception of the distorted signal Y 78, the 
receiver is thus enabled to calculate the channel dependent 
matrix O. The dependency between O and the signals X 77 
65 and Y 78 is expressed in (25) via the matrices R^, the 
autocorrelation matrix of the undistorted signal X 77, Ryy, 
the autocorrelation matrix of the distorted signal Y 78, and 
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Ryx> the crosscorrelation matrix of the signals X 77 and Y 
78. Notice however that, via Rj^ O also depends on the 
delay parameter D. This implies that for each delay value, a 
new channel dependent matrix O must be constructed, each 
of which yields a new Wopt and Bopt. 5 

Simplified models and assumptions may permits the sim- 
plification of the computations necessary to construct the 
channel dependent matrix O. 

In a first simplified implementation of the present 
invention, the undistorted signal X 77 has a flat power 10 
spectral density as shown in FIG. 3A step 114, while noise 
N 76 on the transmission line is negligible as shown in FIG. 
3A step 118. Tims, the correlation matrices RjwCk) and 
R^k) are as follows: 

15 

Therein, 

S JP =E(X / 2 ), 20 
8(k)-l for k-0, and 
6(k)-0 for k*0. 

The autocorrelation matrix R^ therefore becomes equal to: w(i) 

wherein 1 represents the identity matrix. Substituting these 
expressions for R^swOO and Rj^(k) in the above expressions 
(15) for RixOO and Ryy(k), and substituting all the corre- 
. lation matrices in the channel dependent matrix O (25), it is 
found that: 30 



Although the above described embodiments exhibit a very 
low mean square error value and accordingly a negligible 
intersymbol interference outside the guardband, the appli- 
cability thereof in ADSL is doubtful. As already indicated in 
the introduction, a concentration of the energy in the unused 
frequency band must be avoided. The main reason for this 
can best be explained by a frequency domain approach. 
Suppose that H(i) represents the channel transfer function in 
the frequency domain, and E(i) represents an error, equiva- 
lent to E, built up in the frequency domain. Striving towards 
a zero error 

£(0-o (37) 
implies that 

W(l).Y(C)-B(i)X(i) (38) 

wherein W(i), Y (i), B(i) and X (i) are frequency domain 
representations of the equalizer, the distorted signal, the 
target response, and the signal respectively. Equation (38) is 
equivalent to: 



25 and thus to: 

since the channel response H(i) is equal to: 



O - S^fheightf/ • 



(32) 



(39) 



(40) 



(41) 



with: 



Ht*t= 



and: 



in (33): 



^(0) ... r k (T~l) 
rA(r-l) ... r„(0) 

ho ... hx>u-i 



(33) 



If W(i), H(i) and B(i) are split up into unused band and used 
3S band matrices, the above expression (40), because of the 
FDM (Frequency Division Multiplexing) filters used in 
ADSL is rewritten as: 



(34) 



-\ 0 V 



Bowsed 



This results in three conditions: 



-o 



(42) 



(43) 
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M-l 



(35) 



When measuring in an initial step the line impulse response 
vector H, the matrix O can be constructed. 

In a second simplified implementation of the present 
method, a transmission line is also noiseless as in FIG. 3A 
step 118, but the power spectral density of the transmitted 
signal X 77 is no longer flat. The given power spectral 
density of the undistorted signal X 77 must has to be taken 
into account which leads to a new expression for the matrix 
O. 

In another implementation of the present method, noise N 
on the transmission line is not negligible but the power 
spectral density of this noise N is known via additional 
measurements previously executed as in FIG. 3 A step 122. 
In this implementation, a new channel dependent matrix O 
is to be constructed wherein the noise influence is expressed 
via a term R„„(k) which is added to the expression in 
formula (17). R yy (k) is now given by: 

R )V (k>Rj.ky(h k *h^R m (k) (36) 



50 



60 



1W(0^ is unconstrained since B(i) M 
unit energy constraint expressed in the frequency domain 

Although the energy unity norm avoids convergence 
towards the full zero solution, one notices that the unity 
norm does not guarantee avoidance of a full zero solution in 
the useful frequency band. Therefore, the equalizer coeffi- 
cients W must be constrained in order to avoid a concen- 
tration the time domains energy in the used bandwidth. 

In the time domain approach of the above paragraphs, the 
time domains energy can be concentrated in the used band- 
width by adding an additional degree of freedom. The higher 
the energy in the unused frequency band EDC, the less 
optimal the set of equalizer parameters W becomes. This is 
expressed by adding the energy in the unused frequency 
band EDC as an additional contribution to the error E. In this 
way, one obtains the constrained mean square error CE as 
shown in FIG. 4 step 128: 



CE~E+\lE d 



(44) 



Herein, E represents the error contribution due to the dif- 
ference between the equalized impulse response EIR labeled 
Zl 80 and the target impulse response TIR labeled 72 84 on 



07/28/2004, EAST Version: 1.4.1 



with 



-ri-nb -rixb-cr-i) 



5,870,432 

15 16 

FIG. 4, and E^, 238 represents the energy in the unused Still another remark is that intelligent implementations of 
frequency band or don't care band. The suppression factor p. the present should take into account that the calculated 
is a factor that determines how much the energy in this equalizer coefficients will be rounded. Therefore, it is better 
unused frequency band must be suppressed. Writing in to reduce the dynamic range and thus to have higher gains 
matrix notation, one obtains: 5 in the used bandwidth. As a consequence, when different 

criteria provide equivalent solutions for Wopt and Bopt, 
E DC = w T • f£ c • Foe ■ W « w r ■ Q ■ w (45) precision requirements will be taken into account to select 

one of the solutions. By selecting the solution with the least 
dynamic range, the precision performance will be increased. 
(46) io As already noticed above, the channel dependent matrix 
0 depends on the delay parameter D. Thus, the delay value 
which minimizes the intersymbol interference has to be 
found to select the best Wopt. The complete theoretical 
search range for D would be from zero to M. However, 
15 without significant loss of performance, the search range for 
D can be reduced. To reduce the processing time, one could 
Herein, a and b are indexes belonging to the unused or search for example for D values a UmiteA ^tcwi 

"don't care" band. If (13) and (45) are integrated in (44), the arQund Dpea ^ Dpeak being the ^ number ^ G&n 0 
following result is found: and ^ of ^ peak valuc of me transmission line impulse 

E [xr w r .(R YY ^LQ)+B T Ji xx £-2. B r R rx T .W (47) 20 response H as shown in FIG. 3 A step 112 and FIG. 4 block 

112 

It is noticed that the matrix in (13) has been replaced by Furthermore, it is noted that, to reduce the computation 
R^Q. As a consequence, the eigenvector solution repre- u ^ method ^ fee ^ tQ cakukte ^ Qnl 

senting the optimal set of equalizer parameters is found by ei tor w t which has to be foimd m the present 
solving the eigenvalue-eigenvector problem for a new chan- %$ m ^ Q& ^ yefy ^ and mc±Qd is weU 

nel dependent matrix O': but extracts onlv tne eigenvector belonging to the largest 

0'**(Ryy+vQ)-XYx-Rxx~ 1 Yx' W eigenvalue of a matrix. The above eigenvalue-eigenvector 

^jj.-i. - l t. *u * problem (30) therefore is rewritten as follows: 

The additional term p.Q herein can be shown to be equiva- y v J 

lent to the insertion of a constant virtual noise power spectral o.w^KW^O'Kw^.w (50) 

density component into the unused frequency band. The 

term "virtual" is used here to indicate that this is not real _ . r . ■ « 1^ 

measured noise but mathematically injected ooise. ™ eigenvector of matrix O according to. eigenvalue \ 

It is remarked here'that although the described embodi- » *° ™ «f n !f ctor . of ma , tnx ° , aCCOr ^ 0 * 6 

ment in the above paragraphs is a time domain realization, Rvalue X . The eigenvector Wopt according to the 

*i_ 4- *• • i 1- ui ■ *™ 35 minimum eigenvalue Xmm of O is thus found as the eigen- 

the present invention is also applicable in frequency domain 5. . . . _j f . 

r . e iL .( j tt_ vector according to the maximum eigenvalue Amin or the 

approaches of the mean square error method. The error e* r: <^ lulll & w & 

can be rewritten in terms of circular convolutions instead of ma Jf lx . ' . , „ 1 . .. . - „ 

linear convolutions when using the fast fourier transform. In , ™« * e ™ th ° d ^^^r^ZS^rZl 

this way, a frequency domain error criterion is obtained. e, f e ? ^ as shown in FIG. 3Aand FIG. 3B includes the steps 

Such a frequency domain approach will not be described in 40 ° f definug the Ungrt. L of the target impulse response TIR 

further detail in the present application. However, via the ste P m > and defining a P"^*^.^ D 

above described time domain example it is clear to a person 1M ; fo P^. 1 f? 5 ' ,he dlslorted ^ 

skilled in the art how to adapt e.g. the frequency domain » ^l^t k £ t T^f ? H 

approach of J.S. Chow, J.M. Ciofli and JA.C. Bingham in s,e P 120 could be a d, S ltal waveshaping circuit or if deemed 

their article 'Equalizer Training Algorithms for Multicarrier 45 unnecessary for some reason the determining means for Y jn 

Modulation Systems', Proc. IEEE, p. 761-765, 1993 to be ste P I 20 <*> M be nothm S m ° re * an » e , le< f ™ d n ° d ^ A . 

in accordance with the present invention. Indeed, such a predefined sequent representing toe undistorted signd X of 

frequency domain approach can also be adjustedto theFDM S,6 P l 28 ' dis ^ d * e ^ 

M . « j-r ■ lL r j • impulse response TIR and the length of the delay D are used 

requirements by modifying the frequency domain mean ""F«»^ ^ ^ . f 

c/^ • f * ^Jr„;„ 50 to generate the channel dependent matrix O m step 130. The 

square error E(i) in a frequency domain constrained mean . & 1J ^* 1 ,i nt . 

s uare error CEftV channel dependent matrix O is then used to calculate the 

square error c\i). minimum eigenvalue "k MIN of O in step 132 and the asso- 

CE(0 - E{!) + fi{i) • |W(0 2 l ( 49 ) ciated eigenvector of O in step 134. The eigenvector of step 

132 is equivalent to the equalized time domain vector W as 

with 55 determined in step 136. 

A device to calculate a set of equalizer parameters 66 is 

,,(,-)- / ^ if/belon 8 s to ^ WBUBd band shown in FIG. 4 and includes a determining means which 

\ o if i belongs to the used band receives a distorted signal Y 78. The delay D is calculated by 

a means of determining D 104 discussed above and length 
It is further remarked that, for implementation issues, it is 60 L is provided by the definition means of determining TIR 
better to avoid injecting noise close to the FDM filter length (L) 102. The channel dependent matrix (O) is gen- 
boundaries in order to avoid too much attenuation of the first erated by the generation means of channel dependent matrix 
tones of the used band. It is therefore recommended to leave (O) 130 using a predefined sequence representing the undis- 
a buffer zone without virtual injected noise between the torted signal (X) of block 128, the distorted signal (Y), the 
unused and used frequency bands. As an alternative, one 65 delay (D), and the length of the target impulse response 
could use a shaped virtual injected noise instead of a (TIR) (L). The calculation means of minimum eigenvalue 
constant virtual noise. ^mjn °f O anc ^ tne calculation means of eigenvector of 0 134 
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determine the corresponding eigenvector of O. The genera- 
tion means of T equalized parameters W 136 equates the 
eigenvector of O to the equalized time domain vector W 
which is then output from the device. 

As shown in FIG. 5 the time domain equalizer TEQ 16 5 
includes a device to calculate a set of equalizer parameters 
W 66 and an equalizer 90. The distorted signal 34 Y 78 is 
applied to the TEQ and a time domain equalized signal 100 
is output as shown in FIG. 1 and FIG. 5. 

While the principles of the invention have been described 10 
above in connection with specific apparatus, it is to be 
clearly understood that this description is made only by way 
of example and not as a limitation on the scope of the 
invention. 

I claim: 15 

1. A method for calculating a set of T equalizer parameters 
(W), T being a positive integer number, for an equalizer 
included in a receiver to which a transmitter is coupled via 
a transmission line and which is able to receive a distorted 
signal (Y) distorted due to transmission over said transmis- 
sion line of a signal (X), and which is adapted to equalize a 20 
transmission line impulse response (LIR) of said transmis- 
sion line and to thereby generate an equalized impulse 
response (EIR) which approximates a target impulse 
response (TIR) with a predetermined target impulse 
response length (L) delayed over a predetermined delay (D), 25 
said method comprising the steps of: 

a. transmitting said signal (X) and determining said dis- 
torted signal (Y) in said receiver; 

b. defining a first integer number of sample periods which 
represents said target impulse response length (L); 30 

c. defining a second integer number of sample periods 
which represents said predetermined delay (D); 

d. generating a channel dependent matrix (O) from said 
signal (X), said distorted signal (Y), said target impulse 
response length (L) and said predetermined delay (D), 35 
said channel dependent matrix (0) being related to an 
error function that has to be minimized; 

e. calculating a minimum eigenvalue (k MlN ) of said 
channel dependent matrix (O); 

f. calculating an eigenvector of said channel dependent 40 
matrix (0), said eigenvector corresponding to said 
minimum eigenvalue Q^ MIN ); 

g. generating said set of T equalizer parameters (W) from 
said eigenvector, characterized in that said error func- 
tion comprises a first contribution, representing a dif- 45 
ference between said equalized impulse response (EIR) 
and said target impulse response (TIR), and a second 
contribution representing an amount of energy (E^) 
transmitted in unused frequency bands, and further 
characterized in that, to generate said channel depen- 50 
dent matrix (O) said transmission line is supposed to be 
noiseless. 

2. A method, for calculating a set of T equalizer param- 
eters (W), T being a positive integer number, for an equalizer 
included in a receiver to which a transmitter is coupled via 
a transmission line and which is able to receive a distorted 
signal (Y) distorted due to transmission over said transmis- 
sion line of a signal (X), and which is adapted to equalize a 
transmission line impulse response (LIR) of said transmis- 
sion line and to thereby generate an equalized impulse 
response (EIR) which approximates a target impulse <so 
response (TIR) with a predetermined target impulse 
response length (L) delayed over a predetermined delay (D), 
said method comprising the steps of: 

a. transmitting said signal (X) and determining said dis- 
torted signal (Y) in said receiver; 65 

b. defining a first integer number of sample periods which 
represents said target impulse response length (L); 
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c. defining a second integer number of sample periods 
which represents said predetermined delay (D); 

d. generating a channel dependent matrix (O) from said 
signal (X), said distorted signal (Y), said target impulse 
response length (L) and said predetermined delay (D), 
said channel dependent matrix (O) being related to an 
error function that has to be minimized; 

e. calculating a minimum eigenvalue (K MJ ^) of said 
channel dependent matrix (O); 

f. calculating an eigenvector of said channel dependent 
matrix (O), said eigenvector corresponding to said 
minimum eigenvalue (^/at); 

g. generating said set of T equalizer parameters (W) from 
said eigenvector, characterized in that said error func- 
tion comprises a first contribution, representing a dif- 
ference between said equalized impulse response (EIR) 
and said target impulse response (TIR), and a second 
contribution representing an amount of energy (E DC ) 
transmitted in unused frequency bands, and further 
characterized in that said channel dependent matrix (O) 
comprises a first contribution, corresponding to said 
first contribution of said error function, said first con- 
tribution of said channel dependent matrix (O) being 
equal to ^ Y ^^yx ^xx~ 1 '^yx > wherein Ry^ represents 
an average autocorrelation matrix of a sequence of 
samples of said distorted signal (Y), said sequence 
having a length equal to T, wherein R^y represents an 
average autocorrelation matrix of a sequence of 
samples of said signal (X), said sequence having a 
length equal to said target impulse response length, 
wherein R^ represents an average cross correlation 
matrix of said sequence of samples of said signal (X) 
and said sequence of samples of said distorted signal 
(Y), and wherein R^v 1 represents an inverse matrix of 
Rxy and Ryx T represents a transposed matrix of Ryx- 

3. A method according to claim 2, characterized in that 
Ryy(k), which is an element in said autocorrelation matrix 
Ryy whose row index and column index are different from 
each other over an integer value k, is equal to R A:y (k) rh 
(k)+R AW (k), wherein R^(k) represents an element in said 
autocorrelation matrix R^ whose row index and column 
index are different over said integer value k, wherein rh(k) 
represents a convolution of said transmission fine impulse 
response (LIR) and said transmission line impulse response 
(LIR) delayed over k samples, and wherein R^v^k) repre- 
sents an element in a matrix R jVJV , whose row index and 
column index are different over said integer value k, said 
matrix R^ being an average autocorrelation matrix of noise 
(N) measured on said transmission line. 

4. A method for calculating a set of T equalizer parameters 
(W), T being a positive integer number, for an equalizer 
included in a receiver to which a transmitter is coupled via 
a transmission line and which is able to receive a distorted 
signal (Y) distorted due to traasmission over said transmis- 
sion fine of a signal (X), and which is adapted to equalize a 
transmission line impulse response (LIR) of said transmis- 
sion line and to thereby generate an equalized impulse 
response (EIR) which approximates a target impulse 
response (TIR) with a predetermined target impulse 
response length (L) delayed over a predetermined delay (D), 
said method comprising the steps of: 

a. transmitting said signal (X) and determining said dis- 
torted signal (Y) in said receiver; 

b. defining a first integer number of sample periods which 
represents said target impulse response length (L); 

c. defining a second integer number of sample periods 
which represents said predetermined delay (D); 

d. generating a channel dependent matrix (O) from said 
signal (X), said distorted signal (Y), said target impulse 
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response length (L) and said predetermined delay (D), 
said channel dependent matrix (O) being related to an 
error function that has to be minimized; 

e. calculating a minimum eigenvalue °f sa id 
channel dependent matrix (O); 5 

f. calculating an eigenvector of said channel dependent 
matrix (0), said eigenvector corresponding to said 
minimum eigenvalue (X^fflv); 

g. generating said set of T equalizer parameters (W) from 
said eigenvector, characterized in that said error func- 10 
tion comprises a first contribution, representing a dif- 
ference between said equalized impulse response (EIR) 
and said target impulse response (TIR), and a second 
contribution representing an amount of energy (E^) 
transmitted in unused frequency bands, and further 15 
characterized in that, to generate said channel depen- 
dent matrix (O), said signal (X) applied to said trans- 
mission line by said transmitter is supposed to have a 
flat power spectral density. 

5. A method for calculating a set of T equalizer parameters 
(W), T being a positive integer number, for an equalizer 
included in a receiver to which a transmitter is coupled via 
a transmission line and which is able to receive a distorted 
signal (Y) distorted due to transmission over said transmis- 
sion line of a signal (X), and which is adapted to equalize a 
transmission line impulse response (LIR) of said transmis- 25 
sion line and to thereby generate an equalized impulse 
response (EIR) which approximates a target impulse 
response (TIR) with a predetermined target impulse 
response length (L) delayed over a predetermined delay (D), 
said method comprising the steps of: 30 

a. transmitting said signal (X) and determining said dis- 
torted signal (Y) in said receiver; 

b. defining a first integer number of sample periods which 
represents said target impulse response length (L); 

c. defining a second integer number of sample periods 35 
which represents said predetermined delay (D); 

d. generating a channel dependent matrix (0) from said 
signal (X), said distorted signal (Y), said target impulse 
response length (L) and said predetermined delay (D), 
said channel dependent matrix (O) being related to an 40 
error function that has to be minimized; 

e. calculating a minimum eigenvalue of said 
channel dependent matrix (O); 

f. calculating an eigenvector of said channel dependent 
matrix (0), said eigenvector corresponding to said 45 
minimum eigenvalue Q* MIN ); 

g. generating said set of T equalizer parameters (W) from 
said eigenvector, characterized in that said error func- 
tion comprises a first contribution, representing a dif- 
ference between said equalized impulse response (EIR) 50 
and said target impulse response (TIR), and a second 
contribution representing an amount of energy (E^) 
transmitted in unused frequency bands, and further 
characterized in that, to generate said channel depen- 
dent matrix (O), said transmission line is not supposed 55 
to be noiseless and that said method further includes a 
step of measuring noise (N) on said transmission line 
before said channel dependent matrix (0) is generated. 

6. A method for calculating a set of T equalizer parameters 
(W), T being a positive integer number, for an equalizer 
included in a receiver to which a transmitter is coupled via 60 
a transmission line and which is able to receive a distorted 
signal (Y) distorted due to transmission over said transmis- 
sion line of a signal (X), and which is adapted to equalize a 
transmission line impulse response (LIR) of said transmis- 
sion line and to thereby generate an equalized impulse 65 
response (EIR) which approximates a target impulse 
response (TIR) with a predetermined target impulse 
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response length (L) delayed over a predetermined delay (D), 
said method comprising the steps of: 

a. transmitting said signal (X) and determining said dis- 
torted signal (Y) in said receiver; 

b. defining a first integer number of sample periods which 
represents said target impulse response length (L); 

c. defining a second integer number of sample periods 
which represents said predetermined delay (D); 

d. generating a channel dependent matrix (O) from said 
signal (X), said distorted signal (Y), said target impulse 
response length (L) and said predetermined delay (D), 
said channel dependent matrix (0) being related to an 
error function that has to be minimized; 

e. calculating a minimum eigenvalue (h M jx) of said 
channel dependent matrix (O); 

f. calculating an eigenvector of said channel dependent 
matrix (O), said eigenvector corresponding to said 
minimum eigenvalue Qi MZN ); 

g. generating said set of T equalizer parameters (W) from 
said eigenvector, characterized in that said error func- 
tion comprises a first contribution, representing a dif- 
ference between said equalized impulse response (EIR) 
and said target impulse response (TIR), and a second 
contribution representing an amount of energy (E^) 
transmitted in unused frequency bands, and further 
characterized in that said predetermined delay (D) 
depends upon amplitude and index number of a sample 
of said transmission line impulse response (LIR) with 
maximum amplitude. 

7. A method for calculating a set of T equalizer parameters 
(W), T being a positive integer number, for an equalizer 
included in a receiver to which a transmitter is coupled via 
a transmission line and which is able to receive a distorted 
signal (Y) distorted due to transmission over said transmis- 
sion line of a signal (X), and which is adapted to equalize a 
transmission line impulse response (LIR) of said transmis- 
sion line and to thereby generate an equalized impulse 
response (EIR) which approximates a target impulse 
response (TIR) with a predetermined target impulse 
response length (L) delayed over a predetermined delay (D), 
said method comprising the steps of: 

a. transmitting said signal (X) and determining said dis- 
torted signal (Y) in said receiver; 

b. defining a first integer number of sample periods which 
represents said target impulse response length (L); 

c. defining a second integer number of sample periods 
which represents said predetermined delay (D); 

d. generating a channel dependent matrix (O) from said 
signal (X), said distorted signal (Y), said target impulse 
response length (L) and said predetermined delay (D), 
said channel dependent matrix (O) being related to an 
error function that has to be minimized; 

e. calculating a minimum eigenvalue (k MIN ) of said 
channel dependent matrix (O); 

f. calculating an eigenvector of said channel dependent 
matrix (O), said eigenvector corresponding to said 
minimum eigenvalue (Xa/jw); 

g. generating said set of T equalizer parameters (W) from 
said eigenvector, characterized in that said error func- 
tion comprises a first contribution, representing a dif- 
ference between said equalized impulse response (EIR) 
and said target impulse response (TIR), and a second 
contribution representing an amount of energy (E^) 
transmitted in unused frequency bands, and further 
characterized in that before executing steps f and g, 
steps c, d, and e are repeatedly executed, said second 
integer number of sample periods which represents said 
predetermined delay (D) being different in successive 
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executions of step c but laying within an interval of 
acceptable delay values, said channel dependent matrix 
(0) generated in successive executions of step d as a 
result thereof being different in such a way that a 
plurality of minimum eigenvalues similar to said mini- 
mum eigenvalue (^ MI ^) is calculated in successive 
executions of said step e, and that to perform steps f and 
g, only a said channel dependent matrix (0) with 
smallest said minimum eigenvalue (K MJN ) is selected. 

8. A device for calculating a set of T equalizer parameters 
(W), T being a positive integer number, for an equalizer 
included in a receiver to which a transmitter is coupled via 
a transmission line and which is able to receive a distorted 
signal (Y) distorted due to transmission over said transmis- 
sion line of a signal (X), and which is adapted to equalize a 
transmission line impulse response (LIR) of said transmis- 
sion line and to thereby generate an equalized impulse 
response (EIR) which approximates a target impulse 
response (TIR) with a predetermined target impulse 
response length, when delayed over a predetermined delay 
(D), said device including: 

a. a first definition means to define a first integer number 
of sample periods which represents said target impulse 
response length (L); 

b. a second definition means to define a second integer 
number of sample periods which represents said pre- 
determined delay (D); 

c. a first generation means to generate a channel depen- 
dent matrix (O) from signal (X), said distorted signal 
(Y), said target impulse response length (L) and said 
predetermined delay (D), said channel dependent 
matrix (O) being related to an error function which has 
to be minimized; 

d. a first calculation means to calculate a minimal eigen- 
value (k MIN ) of said channel dependent matrix (O); 

e. iteration control means adapted to control said second 
definition means to define a different said second 
integer number of sample periods laying within an 
interval of acceptable delay values in successive 
iterations, adapted to control said first generation 
means to generate a different said channel dependent 
matrix (O) in successive iterations as a result thereof, 
and adapted to control said first calculation means to 
calculate a plurality of minimum eigenvalues similar to 
said minimum eigenvalue Q^ MJJ ^) in successive itera- 
tions; 

f. selection means to select only a said channel dependent 
matrix (O) with smallest said minimum eigenvalue 
Q^min) amongst said plurality of minimum eigenvalues; 

g. a second calculation means to calculate an eigenvector 
of said channel dependent matrix (O), said eigenvector 
corresponding to said minimal eigenvalue (X-vr/zv)? 

h. a second generation means to generate said set of T 
equalizer parameters (W) from said eigenvector, char- 
acterized in that said error function comprises a first 
contribution, representing a difference between said 
equalized impulse response (EIR) and said target 
impulse response (TIR), and a second contribution 
representing an amount of energy (E^) transmitted in 
unused frequency bands. 

9. A demodulator (DEMOD) including between an input 
(DI) and an output (DO) thereof, the cascade connection of 
an analog to digital converter (ADC), a time domain equal- 
izer (TEQ), a serial to parallel converter (SPC) a fast fourier 
transform processing unit (FFT), a firequency domain equal- 
izer (FEQ) and a demapper (DM), said analog to digital 
converter (ADC) being provided to convert a disturbed input 
signal applied to said demodulator input (Dl) into a time 
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domain digital signal, said time domain equalizer (TEQ) 
being provided to digitally filter said time domain digital 
signal and to thereby generate a time domain equalized 
signal, said serial to parallel converter (SPC) being provided 
to transform said time domain equalized signal into a time 
domain parallel sequence, said fast fourier transform pro- 
cessing unit (FFT) being provided to process said time 
domain parallel sequence and to thereby generate a fre- 
quency domain parallel sequence, said firequency domain 
equalizer being provided to digitally filter said frequency 
domain parallel sequence, and, said demapper (DM) being 
provided to demodulate said firequency domain parallel 
sequence and to generate therefrom an output signal applied 
to said demodulator output (DO), said time domain equal- 
izer (TEQ) including a device for calculating a set of T 
equalizer parameters (W), T being a positive integer number, 
for said time domain equalizer (TEQ), said demodulator 
(DEMOD) being coupled to a transmitter via a transmission 
line and being able to receive a distorted signal (Y) distorted 
due to transmission over said transmission line of a signal 
(X), and which is adapted to equalize a transmission fine 
impulse response (LIR) of said transmission line and to 
thereby generate an equalized impulse response (EIR) which 
approximates a target impulse response (TIR) with a pre- 
determined target impulse response length (L), when 
delayed over a predetermined delay (D), said device includ- 
ing: 

a. a first definition means to define a first integer number 
of sample periods which represents said target impulse 
response length (L); 

b. a second definition means to define a second integer 
number of sample periods which represents said pre- 
determined delay (D); 

c. a first generation means to generate a channel depen- 
dent matrix (O) from said signal (X), said distorted 
signal (Y), said target impulse response length and said 
predetermined delay (D), said channel dependent 
matrix (O) being related to an error function which has 
to be minimized; 

d. a first calculation means to calculate a minimal eigen- 
value (k MIN ) of said channel dependent matrix (O); 

e. iteration control means adapted to control said second 
definition means to define a different said second 
integer number of sample periods laying within an 
interval of acceptable delay values in successive 
iterations, adapted to control said first generation 
means to generate a different said channel dependent 
matrix (O) in successive iterations as a result thereof, 
and adapted to control said first calculation means to 
calculate a plurality of minimum eigenvalues similar to 
said minimum eigenvalue (X MW ) in successive itera- 
tions; 

f . selection means to select only a said channel dependent 
matrix (O) with smallest said minimum eigenvalue 
Q*mjn) amongst said plurality of minimum eigenvalues; 

g. a second calculation means to calculate an eigenvector 
of said channel dependent matrix (0), said eigenvector 
corresponding to said minimal eigenvalue (X^^); 

h. a second generation means to generate said set of T 
equalizer parameters (W) from said eigenvector, char- 
acterized in that said error function comprises a first 
contribution, representing a difference between said 
equalized impulse response (EIR) and said target 
impulse response (TIR), and a second contribution 
representing an amount of energy (E DC ) transmitted in 
unused frequency bands. 
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